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1. INTRODUCTION 


The number of physical problems which can be analytically solved is indeed limited. In 
fact, due to complexity of real systems such strict solutions are by their idealisation only 
approximations of reality. Thus finding effective and universal approximate methods is of 
vital importance. In this paper we propose an iterative algorithm of solving the integral * 
- eigenvalue equation for a Hamilton function in the scheme of deformation quantization. 
Our construction is based on the Wentzel - Kramers - Brillouin approximation (WKB) 
developed in the Hilbert space formulation of quantum mechanics. 

The Schroedinger equation can be solved in an exact way when the form of the potential 
is relatively simple. For other cases several approximated methods have been developed 
during the years. One of such techniques is the WKB approximation mentioned in the 
previous paragraph and sometimes called the quasi-classical approximation. This 
method was introduced in 1926 and is suitable for quantum systems where the potential 
changes slowly in comparison to the de Broglie wavelength. 

The WKB algorithm works well for the position representation. The main difficulty in its 
application arises from the fact that the domain of a solution of the Schrodinger equation 
must be divided in spatially separable regions. On each region one gets a piece of the 
wave function and the problem is to put these pieces together to obtain a global solution 
over all the regions. This obstacle is overcome with special connection formulas. 

On the other hand, the phase space quantum description is an alternative approach to 
the Hilbert space formalism of quantum mechanics. This scheme provides supplementary 
valuable information in addition to the usual quantum formulation that is carried out in 
just one representation (position or momentum). In this framework the coordinates and 
momenta are considered simultaneously, which gives a natural extension of the Hamilto¬ 
nian construction to describe quantum systems. A complete review of this topic can be 
found in [6HH] - 

In the phase space formulation of quantum mechanics the information about the system 
is obtained through a Wigner function, which plays a similar role as the wave function in 
usual quantum mechanics. However, by means of this function it is possible to study the 
classical limit in a more transparent way. The Wigner function has been employed widely 
in quantum optics, condensed matter, nuclear and particle physics, etc. A comprehensive 
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guide to these applications can be found in [7] and the references cited therein. 

A general way to obtain the phase space quantum description of a system is through the 
deformation quantization formalism. Its main advantage is that it can be used to treat 
systems with arbitrary phase spaces. To obtain a Wigner function for an arbitrary system 
it is necessary to find solutions of the so called * - eigenvalue equation. However this 
integral equation is in general difficult to be solved and, as with other physical theories, 
certain approximation methods are needed to obtain a solution. 

Our intention is to propose an adaptation of the WKB method in order to obtain ap¬ 
proximate solutions of a * - eigenvalue equation for a Hamilton function in the phase 
space description of quantum physics. In this approach, due to the fact that the positions 
and momenta coordinates are used on an equal footing, a straightforward recalculating of 
the quasi - classical formulas fails and some additional considerations are required. For 
example, the Weyl correspondence connecting wave functions and their respective Wigner 
functions is nonlocal. Thus we are faced to the problems of representing a product and a 
sum of wave functions in terms of the respective Wigner functions. Solutions of these two 
problems are analysed in this article. We also focus attention on an important question of 
representing an interference of wave functions as a contribution to the Wigner function. 
The considerations presented in this work can be related to the problem of a semiclassical 
limit of a Wigner function ®-H- However, we do not have to handle with singularities 
of Wigner functions, because difficulties caused by the term | are eliminated at the level 
of wave functions. It is worth to mention that in another context an application of the 
WKB expansion in formal deformation quantization has been studied in [HI ' 15] . 

Our paper is organised as follows. In Sec. 2 the energy eigenvalue problem in the defor¬ 
mation quantization formalism is briefly described as well as the most important aspects 
of this quantization approach. Then, in Sec. 3 we present a review of the main elements of 
the WKB construction. Next, in Sec. 4, we transform the quasi-classical approximation 
for wave functions into an approximation for the respective Wigner functions. Specifically, 
we propose formulas to represent a product of wave functions and a sum of wave func¬ 
tions with separate supports. We introduce a Wigner function representing interference 
of states and analyse some of its properties. Two examples of this WKB construction 
are presented in Sec. 5. General considerations are illustrated by calculations done for 
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a one dimensional (1-D) harmonic oscillator and for unbound states in the Poeschl - 
Teller potential. There are also two Appendices. The first one contains a derivation 
of the Wigner function for states represented by tempered distributions. In the second 
appendix the Wigner functions of unbound states in the Poeschl - Teller potential have 
been calculated. Finally, we give our concluding remarks. 


2. THE ENERGY EIGENSTATES PROBLEM IN DEFORMATION QUANTI¬ 
ZATION 


The phase space quantum description is an alternative approach to the Hilbert space 
formalism of quantum mechanics. In this formalism the corresponding space of states is a 
symplectic manifold and observables are represented by smooth real functions. Thus the 
coordinates and momenta are considered simultaneously, which gives a natural extension 
of the Hamiltonian construction to describe quantum systems. The idea is based on the 
fundamental observation that quantum physics is a deformed version of classical theory 
(see [IB]) where the role of deformation parameter is played by the Planck constant h. 

A general way to obtain the phase space quantum description of a system is through 
the deformation quantization formalism. Its main advantage is that it can be used to 
treat systems with arbitrary phase spaces although, since the multiplication of functions 
is replaced by a new product called the * - product, equations appearing in the formalism 
are in general difficult to solve. 

Nevertheless, the deformation formulation of quantum mechanics liberates us from nu¬ 
merous formal obstacles. First of all we do not need to construct an associated Hilbert 
space of the system. We also avoid quantization of observables and the nontrivial problem 
of defining domains of the constructed operators. In its general version the deformation 
quantization calculus works on nontrivial symplectic spaces whereas other quantization 


procedures may not be established at all. Thus e.g. the * - eigenvalue equation (2.1a) 


and (2.1b) presented below is well defined on an arbitrary symplectic manifold while the 
stationary Schroedinger being its counterpart may be not known. 

Our goal is to obtain an approximated method to find eigenvalues and eigenstates of a 
Hamilton function in the framework of deformation quantization. Thus we start from the 
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* - eigenvalue energy employed in this approach. 

The * - eigenvalue equation for a Hamilton function H (r, p ) has the following form 

H(f,p) * W E (r,p) = E W E (r,p), (2.1a) 

where E denotes an energy eigenvalue and We(t,p) the corresponding Wigner energy 
eigenfunction. Moreover, the additional condition 


{H(r,p),W E (f,p)} M = 0 


(2.1b) 


on the Wigner energy eigenfunction We (r, p ) is imposed 
On the phase space M 6 , as the * - product we use the Moyal product p~7l flTJ] 

1 


A(r,p) * B(f,p) := 




dr'dp'dr"dp"A{f' ,p')B(f", p ") 


x exp 


Y | (f" — r) ■ (p' — p) — (r' — r) ■ (p 1 


P. 


( 2 . 2 ) 


where the dot stands for the scalar product. The above definition of the * - product is 
valid for a wide class of tempered distributions (for details see [20]). This observation is 
important since the Wigner energy eigenfunction We (r, p ) can be a generalised function. 
The sign convention employed in this paper is compatible with the Fedosov works [21, B21. 
The Moyal product is associative but in general non-Abelian. Moreover, it is closed i.e. 

/ A(r,p) * B(r,p)dfdp — / B(r,p) * A(r,p)drdp = / A(rpp) ■ B(r,p)dr dp. 

Jr 6 J r 6 Jr 6 

The Moyal bracket appearing in the condition ( 2.1b[ ) is defined as 

{A(f,p),B(r,p)} M := ^ (A(f,p) * B(r,p) - B{r,p) * A(f,p)^J. (2.3) 

From these definitions it can be deduced that formulas (2.1a) and (2.1b) are integral 
equations and there is no universal way of solving them. Thus our idea is to adapt 
the WKB approximation from the Hilbert space formulation of quantum mechanics to 
obtain approximated Wigner energy eigenfunctions. Such adaptation seems to be possible 
because for systems with phase spaces of the type M 2n a correspondence between the phase 
space description and the Hilbert space description is known. 

Indeed, if we restrict to the case of the phase space M 2 , the phase space counterpart 
A(x,p) of an operator A acting in the Hilbert space L 2 (M) equals 


A(x,p) = W A) = 


“ + OO 


d i ( x -\ 


A 


£ 


x + - ) exp ( —— 




(2.4a) 
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A similar formula in the momentum representation holds 

/*+OO 


A(x,p) = W" 1 (A) = 


d'l (P ~ \ 


A 


77 \ / 777a; 

P+ 2/ eXp X 


(2.4b) 


The mapping W 1 is called the Weyl correspondence (for details see [19]). For every 


operator A such that its representation (x — | 


A 


x+ 1) is a tempered generalised function 


of £, the generalised function A(x,p ) is well defined. Note that the Weyl correspondence 


(2.4a, 2.4b) is nonlocal. 


Applying the Weyl correspondence W _1 to a density operator of an energy eigenstate \iPe) 
with the wave function i(je{x) it can be observed that the respective Wigner function is 
of the form 


W E {x,p) := W 


-l 


\^e)(^e\) = 


2ttH 


2ttH 






h 

(2.5) 


Again the relation between the Wigner function and its respective wave function is nonlo¬ 
cal. The value of We(x,p) at a fixed point (x,p) depends on values of the function i/je(x) 


on the whole configuration space. A strict form of expression (2.5) applied to tempered 
generalised functions can be found in Appendix A. 


3. THE WKB APPROXIMATION IN WAVE QUANTUM MECHANICS [H [2| 


Let us consider the one particle nonrelativistic Schroedinger equation 

: A f) + V(r )^{t, f) = ih 9 ^ r ^ 


(3.1) 


2M T ' ' ■ V 27-V-, w dt ■ 

If the potential V(r) does not depend on time, the stationary wave function ^( t,f) is of 
the form 


r) = exp 


iEt\ 


h 


I V’B(f), 


(3.2) 


where E denotes an eigenvalue of the Hamilton operator and (F) its corresponding 
eigenfunction. Hence ^e(x) obeys the stationary Schroedinger equation 

h 2 


2 M 


At /} E (r) + V(r)i/j E (r) = E^ E {r). 


(3.3) 


Moreover, every solution of (3.3) can be written as a linear combination of two functions 

= exp f jfiif) ) and ip E n(r) = exp f ^cr 77 (f) J (3.4) 
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satisfying separately Eq. (3.3), where oy(r) and ajj(f) denote some complex valued func¬ 


tions. The phases oy(r) and an(r) fulfill the second order nonlinear partial differential 
equation 

1 o ifi 

(3.5) 


:(Vcr(r )) 2 — (r ) = E — V(f), 


2 M y y " 2 M 

for cr(r) = oy(r ) and cr(r ) = <Tn(r). In the classical limit h —* 0 this equation reduces to 
the Hamilton - Jacobi stationary equation 

1 


2 M 


(Vcr(r)) = E — V(f), 


(3.6) 


where the function cr(r) is interpreted as the stationary action and its partial derivatives 


dap Oap aap ^ ^ mQmenta _ 
ox 7 oy 7 oz 


The expression (3.5) is equivalent to the stationary Schroedinger equation (3.3). The 
phases cr/(r), cr//(r) carry the same information about the eigenstate of the Hamilton 
operator as the energy eigenfunction does. However, it is usually more difficult to solve 


(3.5) than Eq. (3.3). Importance of the formula (3.5) lies in an iterative procedure in 
order to derive the function cr(r). 

For simplicity let us consider the 1-D case. Then Eq. (3.5) is of the form 

2 


2 M 


da(x ) 


ih d 2 a(x ) 


= E- V(x). 


(3.7) 


dx J 2 M dx 2 

In certain parts of its domain the solution can be written as a formal power series in the 
Planck constant 


o (x = 


E 

k =0 


°k[X) 


(3.8) 


Thus we receive an iterative system of equations 


2 M 


da 0 (x) 

dx 


dao(x) dai(x) 1 d 2 cr 0 (x) 

dx dx 2 dx 2 


da 0 (x) da 2 (x) 
dx dx 


1 

+ 2 


d(T\{x) 

dx 


1 d 2 cr 1 (x) 

2 dx 2 


= E-V(x), 
= 0, 

= 0, 


(3.9) 


The element cr 0 (x) is the classical stationary action for the system while the next terms 
can be interpreted as quantum corrections. 
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The recurrence (3.9) leads to the following set of conditions 


= ±y/2M(E-V(x)) 
dx 

da n (x) _ 1 „ / d 2 (r n _i(x) da n _i(x) da^x) ^. 

d X da 0 {x) n l d x 2 ’ fa ’ ’ dx ) 

dx v 7 


(3.10) 


By f n (yi, ..., y n ) we denote a polynomial in variables y±,..., y r 


There exist two solutions of Eqs. (3.10) and they differ on the sign at even h power 


elements. Thus the phases crj(x) and ajj[x) from the formula (3.4) are 

°° /1 \ k oo / j- \ k 

ai(?) = ( 7 ) C 7 fc( x ) > ( - ) (-l) fc+ V fc (a:). (3.11) 

k =0 ' 1 7 fc=0 ^ 1 7 

The odd coefficients opfc+i can be selected real and the even elements cr 2 /c k — 0,1,2,... 
can be chosen real for E — V(x) > 0 and imaginary for E — V(x) < 0. This freedom 


of choice results from the fact that the system (3.9) determines only the first derivatives 


of cTjt, k — 0,1, 2,... ffence an arbitrary complex number can be added to each term cq,. 
Thus the energy eigenfunction is normalisable and one can write 


op (a:) = 0- O dd(x) + O’even(^) , CF n (x) = CT odd (x) - CTeven(^), 


(3.12a) 


where 


<7odd(^) := Y. ( j CT2fc+l(ic) , 


k=0 


@( % ) • 


==£ 


OO 2k 


k =0 


CT2k(x). 


(3.12b) 


With respect to the choice: (J 2 k+i( x ) ~ real, cr 2 fc(ar) ~ real or imaginary we conclude that 
crodd (*r) is an imaginary function and cr ev en(^) is a real function for E — V(x) > 0 and an 
imaginary function for E — V(x) < 0. 

The series expansion cannot be applied in neighbourhoods of turning points. The system 


(3.9) consists of infinitely many equations. For obvious purposes we would like to restrict 


ourselves to a finite number of terms in the series (3.8). ffowever, it is possible only if a 


few initial terms make a principal contribution to the sum (3.8). Thus the necessary but 

















not sufficient condition for the expansion (3.8) to be applicable is that it should be taken 
at points x satisfying the inequality 


\x - x 0 | > 2 


h 2 


M 

dV{x) 1 


dx 1 x=xq 


(3.13) 


(for details see [I] [2j). Certainly the turning points never obey (3.13) and their neigh¬ 
bourhoods must be analysed separately. 


One has to remember about another restriction imposed on the expansion (3.8). The 


turning points are places at which the momentum p = ±\/E — V(x) is zero. Hence from 


the first equation of the system (3.9) we see that the expansion (3.8) may not work at 


points at which the classical momentum is small. 

In order to find a complete approximate solution of the stationary Schroedinger equation 


we have to solve first the system of equations (3.9) in all intervals, in which the series 


expansion (3.8) is acceptable. Then we match these approximate solutions. There are 


several techniques of making the separate approximate solutions compatible [4J [5]. The 
one used here applies strict solutions of the stationary Schroedinger equation near turning 
points. 

The algorithm sketched above is widely used and appears in the literature as the quasi - 
classical approximation or the Wentzel - Kramers - Brillouin (WKB) approximation. 

Let us discuss its application to a 1-D potential drawn on FIG. [l] We look for a solu¬ 
tion of the Schroedinger stationary equation using the first order of the quasi-classical 
approximation. 


The domain of Eq. (3.7) has been divided in five regions belonging to three classes. 


1. Classically forbidden regions A L = (—oo,cq>] and A R = [a^oo). In these regions 
the iterative system of equations (|3.9) is applicable. 


2. Neighbourhoods Bl = [&Lz,&i>] and B R = \b R i,b Rr \ of the turning points in which 
is necessary to approximate the potential V (x) by a polynomial and then to find a 


strict solution of a 1-D version of the stationary Schroedinger equation (3.3). 


3. A part C = [q, c r ] of the classically accessible area between the turning points x± 


and x 2 , in which the iterative procedure (3.9) can be used. 
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FIG. 1: A potential V{x) as a function of x. 


In the area Al the energy E is smaller than the potential V(x). From the first equation 


of the system (3.9) we obtain that 

Vo(x) = ±t f P (v)dy , p(y) := V 2M (V(y) - E). 


(3.14) 


'Xl 


The choice of the sign results from the requirement that the spatial probability of detection 
of the particle must be normalisable. Thus the unique physically acceptable first order 


solution of the system (3.9) is the pair of functions 

[•X 

<t 0 (x) = ~i P (y)dy , ci(x) = In 

J Xl 

The respective wave function equals 


n/p(z) 


+ G,, G x e C. 


(3.15) 


^E(i)A L {x) = -y== exp (j I P (y)dy \ , (3.16) 

VPW \ n Jx 1 J 

where D^a l corresponds to a normalising factor. 

For the region Bi the potential V(x) is approximated by a polynomial. Since we consider 
only the first step of the quasi-classical approximation, it is sufficient to assume that the 
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potential is linear 


V{x) « V{x{) - F Xl (x -xi) , F Xl 


dV(x) 


dx 


X=X\ 


In this case the stationary Schroedinger equation reduces to the form 

d 2 ^EB L {x) , 2 MF X1 


+ 


-(x - Xi )i>EB L { x ) = 0. 


(3.17) 


(3.18) 


dx 2 h 2 

The probability of detection the particle diminishes as the coordinate x decreases. Thus 


the unique acceptable solution of (3.18) is 


EB L (x ) — Cb l ® 


^ 2 MF xi y 3 


(xi - x) , 


(3.19) 


where &(y) denotes the Airy function defined by &(y) / 0 °° cos + uyj du and the 

constant Cb l £ C. 

In the intersection area of regions Al and Bl, where X\ — x » 0 and the relationship 


(3.13) holds, an asymptotic expansion of the solution (3.19) can be applied (for details 


see 0)- The asymptotic expansion of ( |3.19[ ) and the function ( |3.16[ ) coincide under the 
condition D^a l — —^(2 MF xl h)^ 6 . 

The region C is classically accessible and the wave function in this area can be chosen 
real so in the first quasi-classical approximation we have 

i>E(i)c(x) = - 7 == sin (If P (y) d y + ^ , P(l/) : = V 2M ( E - V(y)), D (1)c , 5 e M. 

VPW \ n Jx 1 J 

(3.20) 

Inside the intersection of the regions C and Bl, where X\ — x « 0 and the condition 


(3.13) is fulfilled it can be employed the asymptotic expansion of (3.19) 


Therefore we can see that 


S= 7 L + kn, ke Z and D (1)c = (-l) k C BL (2MF xl h) 1/6 . 
If we take the simplest choice k — 0, the wave function is given by 


il>E(i)c(x) = 


D 


(i )c 


7r 


p{v)dy + 4 


'X 1 


(3.21) 


(3.22) 


In a way analogous to the earlier considerations we obtain that in the area Ar the wave 
function equals 


i’E{l)A R (x) = 


D 


k=<*P — z / P (v)dy 

vpM V h Jx 2 


(3.23) 
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In a neighbourhood of the turning point X 2 we have a wave function similar to (3.19) 

dV{x 


/o j\/f r' \ 1/3 

= C7 fla $ ( (*2 - *) ) , F X2 := 


dx 


(3.24) 


X=X 2 


The coincidence condition between the solutions in the areas Ar and Br means that there 
must be D(i)a r = — 1 ^(— 2MF X2 K) 1 ^. 

Finally the wave function in the region C seen from the side of the turning point X 2 is 


i ; E(l)c( x ) — 


D\ 


(1)0 


\Zp(x) 


sin 


1 

h 


P (y)dy + 5'j , p(y) := y/2M(E-V(y)), D[ 1)C ,S ' G R. 

(3.25) 


As before, comparing the approximate solutions in the intersection C fl Br we conclude 
that 

S' = ^ + kn , k E Z and D[ 1)c = (-l) k C BR (-2MF X2 H) 1/6 . (3.26) 

Again, selecting fc = Owe obtain 


' i Pe(i)c( x ) — 


D' 


(1 )c 


71 




' 0)2 


(3.27) 


The functions (3.22) and (3.27) must be equal. Since the integral p (y)dy is positive, 


we arrive to the well known quantization rule 


rx 2 / 1 ' 

P (v)dy = h(n + - ) 7T, neJ\f. 


(3.28) 


The coefficients and fulfill the equality D(i)c = (— l) n D'^ c . 

I 11 the next chapter the elements of the WKB construction presented here will be extended 
to the deformation quantization formalism. 

4. AN ADAPTATION OF THE WKB APPROXIMATION TO DEFORMATION 
QUANTIZATION 


The main idea of an approximation method in the WKB approach to solve the energy 


eigenvalue problem (2.1a, 2.1b) for a Hamilton function H(r,p) is to explore the equa¬ 


tion (3.5) for a phase cr(r). In the 1-D case a relationship between an energy Wigner 


eigenfunction W E (x,p) and the phase a(x) follows from the Weyl correspondence (2.5) 

^ r+00 


W E (x,p ) = 


2 ^ exp U 


° - a (a:+ ^ ) - £p 


(4.1) 
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Notice that the relation between the Wigner function and the phase is global i.e. the 
value of We(x,p) at a fixed point depends on values of the phase a(x) on the whole 
configuration space. 

Unfortunately, as we mentioned in the previous section, it is usually necessary to restrict 


to an approximated solution of Eq. (3.7). In particular, the iterative WKB algorithm of 


solving this problem has been discussed. Thus to propose an effective method of dealing 
with the energy eigenstates problem in deformation quantization we have to transform 
elements of the WKB approximation for wave functions into procedures dedicated to 
Wigner functions. 

In fact four elements must be considered: conditions imposed on eigenvalues of the Hamil¬ 
ton function, representation of the sum phases in terms of phase space functions related 
to these phase, a relationship between a superposition of wave functions and their Wigner 
functions and finally, the connection formulas for partial Wigner functions. 

As it could be seen in Sec. [3j there is no universal rule determining energy eigenvalues 
in the WKB method. For example in the first approximation for a potential like the one 
depicted at FIG. [l] the condition imposed on admissible values of energy is of the form 


(3.28). 


Three other components of the WKB procedure require more detailed analysis. Thus to 
each of them a separate subsection will be devoted. 

4.1. The phase space representation of a product of wave functions 


Let us go back to the formula (3.8), on which the WKB approximation is based. We see 


that in the case when the phase <j(x) is a power series in the Planck constant, the wave 
function has the following form 


^e(x) = W^Ek{x) , 'lpEk(x) ■= exp 


k =0 


fiU 1 atW 


(4.2) 


The functions ip E k{x), k = 0,1,... need not be elements of L 2 (M) but as they are smooth 
and, due to physical requirements, bounded, the product ip Ek (x + |) ip E k {x — |) is a 


tempered generalised function for every k. Since the Weyl correspondence (2.5) is in 


general a Fourier transform of the product of functions ip {x + ip {x — §) (look the 
Appendix [ a|, then the integral d£ ip Ek (x + |) ip E k (x — §) exp (—yf) is well defined. 
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Under the conditions discussed in Sec. [3] the finite product rifc=o ^Ek{x) may be a good 
approximation of an energy eigenstate. This approximation can be realised as an iterative 
procedure, in which the n-th approximation ^E(n){ x ) °f the wave function i])e{x) equals 


lpE{ 0)(x) := 1pEo(x) , fJ}E(n)(x) = ^E(n-l )( x ) ' ^En( x ) , U> 1 . 


(4.3) 


Thus we are interested in finding a formula expressing We(ti)( x ,p) '■= 

W_1 (2S»l^(")X^(")l) by W E{n-i){x,p) := W- 1 ^|^ E(n _ 1) )(^ E(n _ 1) |j and 

WE n (x,P ) := W- 1 . From (2.5) it can be observed that 


f + OO 


iXp 


h 




dpW E (x,p)exp — =^e\x+-)^e[x-- 




X 


This result in fact is true for an arbitrary pure state. Now from (2.5) and (4.3) we obtain 

i ^+00 r +00 r+oo 


W E (n){x,p) = 


27 ih 


dp' 


dp" 


dx W E (n-i) {x, p')W En (x, p")exp 


iX(p' + p" — p) 

h 


r»+QO 


r»+oo 


= / dp'W E (n-l)(x,p')W E n(x,p-p') = / dp"W E {n-l)(x,p-p")W E n{x,p"). (4.4) 

J — oo J —oo 

In this way the Wigner function corresponding to the product of two arbitrary wave 
functions is represented by the convolution of their Wigner functions with respect to the 
momentum. 

Analogously, if we work in the momentum representation with a function of the form 


1>B( 0)(P) := feo0) , fe(,)W = V>E(„-l)(p) ' 1pEn(p) , « > 1, 


(4.5) 


then 


r >+00 


r >+00 


IU B ( n )(x,p) = / dar , W S (n-i)(a: , ,p)WB n (a;-a; , ,p) = 


dx ,, lU E (n-i) (x-x", j?) W Sn (x", p) 

(4.6) 


4.2. The phase space counterpart of a superposition of wave functions 


The main difficulty in the application of the WKB method is that the series expansion 


(3.8) cannot be applied everywhere. Thus in the quasi-classical approximation the total 


wave function is represented by a sum of k spatially separable functions 

k 

V'e(x) - 2j^E Ui 6i( x ) . -oo < < b 2 = a 2 < b 2 = a 3 < ... < b k -i = a k < b k < oo. 

(4.7) 


1=1 
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The function ip E {x) need not be continuous at the points a*, 6./, i — 1,2,... k. Nevertheless, 
every partial function 'ip Eai b l ( x ) satisfies the condition supp^sa^x) C [ai,bi]. Therefore 

ipEa l b l {x) = Y(x - ai)ip Eai b l {x)Y(bi - x) = Y{x - ai)tp E {x)Y ( 6 / - x), 


where Y(x) denotes the Heaviside function. In this way the set of functions 
{'ipEaibi( x )}f=i is orthogonal. Provided the global function iPe(x) is normalised, we obtain 

Ef=l W^EaibXx)]] 2 = 1. 

Dehning the self-adjoint operators 


Pr 


Eaibi 


{fpEaibilfpEaibi) 


| IpEafy ) {i’Eatbi |, l — 1, ■ ■ ■ ,k 


(4.8) 


which are projectors on 1-D closed subspaces of a Hilbert space PL, one can see that their 
product satisfies the relation 


PfEa^PfEa r b r ~ $lrP r Ea r b r - 


It is now of vital importance to obtain a phase space counterpart of the state being the 


superposition of wave functions of the form (4.7). Let us consider a Wigner function 


arising from a wave function / ipEa l b, (x). Applying (2.5) one gets 


W Eaibl (x,p) = 


r'Min.[2(x—ai),2(bi—x)] 


2-Kh 


Max.[2(aj— x),2(x—b[)] 


dpip 


Eaibi 


p 


X + - ) Ip Eaibi ( X - - ) exp (- — 


ipP 


(4.9) 

Therefore the Wigner function WE ai bi(x,p) vanishes outside the set ( ai,bi ) x M. As the 
function ip E aibi(x) itself can be a sum of functions, we observe that every Wigner function 
corresponding to a superposition of wave functions with supports from an interval [ ai,bp\ 
is still limited to the strip a/ < x < 6/. 


Moreover, for a/ < x < the integral (4.9) turns into 


"2 {x-ai) 


dpip 


Ea t bi 


J2(ai-x) 

and for < x <bi we obtain 

f2(bi-x) 


P 


x + 5 ^Eaibi X - - exp —— 


ipP 

h 


' 2(i—6() 


dp ip 


Eaibi 


P 


X + - IpEatbi \X ~ ~ exp-— 


ipp 

h 


The interval of integration always contains 0, is symmetric with respect to the point p = 0 
and its length increases from 0 for x = ai to 2(6/ — a/) for x = Then it decreases to 


15 
















0 for x = bi. As a Wigner function, the function W E aibi{x,p) is continuous with respect 
to x and p and finite at every point of M 2 . If ai and bi are finite, the function WEaibi(x,p) 
as the Fourier transform of a function with compact support, is smooth and its support 
is unbounded. 

Analogously, for every wave function of the form i^Ecid^P) = Y(p — q)?/> e(p)Y( di — p) , q < 
di in the momentum representation its Wigner function WE Cl di(x,p) vanishes for p < q 


and p > di. 


From the expression (2.4a) we can deduce that if an operator A in the position represen¬ 
tation satisfies the condition (x\A\x'} ^ 0 only for a < x,x' < b, then the function 
W -1 (A)(x,p) may be different from 0 only for x contained in the interval (a, b). Moreover, 
the function W _1 (A)(:r,p) is a smooth function with respect to the momentum p. For 
every x G (a, b) and every positive number A > 0 there exists a value of momentum p 
such that \p\ > A and W _1 (A)(x,p) ^ 0. 

An analogous conclusion can be formulated for every operator A such that in the momen¬ 
tum representation (p\A\p'} ^ 0 only for c < p,p' < d. 

It is important to point out that the inverse statement is not true i.e. the fact that a 
function A(x,p) = W _1 (A) differs from 0 only for a < x < b or, respectively c < p < d, 
does not imply that (x\A\x , 'j ^ 0 exclusively for a < x,x' < b or, respectively (p\A\p'^ ^ 0 
exclusively for c < p,p' < d. 

Indeed, let A(x,p) = Y(x + a)Y(a — x)Y(p + c)Y'( c — p) with a, c > 0. Then 


(x\A\x') 


sin 
2 h — 


^ c(x-x') j 
X — X 1 


for —2a < x + x' < 2a and 0 for x, x' not fulfilling these inequalities. 

Finally, for every 1 < l,r < k one finds that 

rW roo 

/ dx dpWEathix.p) = \\lpEatbM 2 

J ai J— oo 

and 

[ dx f dpW Eai bi (a, P)W E arb r ( x , p) = Sir ^2 a h ^ ^ ' 

«/ CL l —OO 

Now we are going to analyse a problem of representing a superposition of wave functions 
on the phase space. Let us consider a two-component linear combination of functions 


Y{x~ai)ipEa l b l {x)Y(bi-x)-\-Y{x-a r )il)Earbr{x)Y(b r -x), -oo < ai <bi < a r <b r < oo. 
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Its Wigner function is then given by 


W E {x,p) — W 1 {^^\i , Ea l b l ){'^Ea l b l \ ) + W 1 ( 7^ \^Ea r b r ) {^Ea r b r | ) + 

+ W_ 1 (^ J ^l' i / ; Ea ! 6 i )(V , J Ea r 6 r | + Ea r b r ) (P’Ea^- ( 4 - 10 ) 

The components W _1 ^^\ipEaibi)(i’Ea l b t \ ) and W _1 (^\ipEa r b r )(i’Ea r b r \^ belong to the 
set of partial Wigner functions analysed before. However the other two terms are essen¬ 
tially new and they correspond to the function 


Wem {x,p) := W Ea r b r I + ^^\'pEa r b r ){'pEa l b l \ ) , (4-11) 

which represents the interference between the partial wave functions ipEaibi(x ) and 

^ Earbr {x). 


The operator appearing in (4.11) and defined as 


lilt |'0-Ea/6;)('0-Ea r & r | + | IpEa r b r ) Eaib/. 


(4.12) 


is self-adjoint although it is not a projector. Its trace vanishes and it has three possible 
eigenvalues A: 



(4.13a) 


A 0 = 0 , its eigenvector is every vector orthogonal to \'ip E a l bi) and |t/iEa r & r )> (4.13b) 


\+ = W^EcHbiW ■ \\i>Ea r br\\ , I T ) ~~7= (JT~j - rrl^Eaibt) + ITT - Tr\lpEa r b r ) ) • (4.13c) 

V 2 WlVEatbiW | \YEa r b r \ \ J 

The interference operator Int exchanges the directions between vectors \ipEaibi) ^ 
| Ea r b r )- Indeed, 

Intl^E^fe,) = I |VtEa;b ; 11 2 \P’Ea r b r ) , hit \^Ea r b r ) = | |'0^7a. T .6 T . 11 2 | 'IpEafo)- 


The function W E int(%,p) defined by the expression (4.11) and representing the interference 
term is determined by the real part of the integral 


W Eint (x,p ) = 2 Re 


r»Min.[2(6j— x),2(x— a r )] 


' Max.[2(aj — x),2(x—b r )] 


d£ip 


Eaibi 


X + | ) Ip Ea r b r 


A / ipp" 

x ~2 exp Hr 


(4.14) 
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It may be different from 0 exclusively for x E ( a; ^ ar , b ’+ br ^ . This interval in general is 
not contained in the sum of intervals ( ai,bi ) U (a r ,b r ). Hence the interference part of the 
Wigner function ( |4. 10 ) can be nonzero only at points with abscissae, at which the wave 
functions i/)E ai bi{x) an d ipEa r b r ( x ) disappear. The function We int,( x , p) is real and it does 
not contribute to the spatial density of probability, because 


r»+oo 


Qint(x) = / dpW E int(x,p) = 0. 


(4.15) 


Therefore 


bi-\-b r 

r+ oo r+oo C~2 — /*+oo 

/ dx dpW Eint (x,p) = / dx dpW Eint (x,p) = 0. 

— OO «/ — oo «/ a ^~f ar «/ —oo 


The following integrals also vanish 


/*+oo /»+oo /»+oo /*+oo 

/ dx dpW E int(x,p)W E a l b l (x,p)= dx dpW Eint (x,p)W E a r b r (x,p) = 0 . 

OO j—OO j —OO j —OO 


Both limits of the integral (4.14) with respect to <i£ are always negative and the length of 


the interval of integration fulfills the inequality 


Max.[2(a/ — x), 2(x — b r )] — Min.[2(6; — x), 2(x — a r )} < Min.[2(6; — a;), 2(6 r — a r )}. 

For any observable A{x) depending only on the position x, the interference Wigner func¬ 
tion W E i n t(x,p) does not influence the mean value of A(x), because 

/ +oo r+oo 

dx dpW Eint (x,p)A(x) = 0. 

oo J — oo 

The existence of the interference Wigner function leads to an apparent paradox. Indeed, 
assume that a wave function, which is not necessarily an energy eigenfunction, is of the 
form 


4>(x) = y(z - ai)ip ai b 1 (x)Y(bi -x) + Y(x - a 2 )^ a2 b 2 ( x )Y ( 6 2 - x), 


(4.16) 


—oo < di < bi < a 2 < 6 2 < oo. 


Thus this function vanishes between points b\ and a 2 - Therefore the density of probability 
of detection of the system in any spatial point from the interval [ 61 , 02 ] is zero. From our 
previous considerations we deduce, that for any observable A(x), = 0 with the 
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support being a strip in the spatial interval [&i, a 2 ], the mean value (^A(x)J> = 0 although 
the Wigner function W(x,p) of this state is different from zero there. 

Even a stronger conclusion can be formulated. Let A(x,p) be a real function defined as 


£) ex P (“If) for h<x<a 2 


A(x,p) = < 


(4.17) 


0 for x < b\ and x > a 2 . 

The function f(x,£) is chosen in such a way that it ensures reality of the function A(x,p) 


and the existence of the integral in the formula (4.17) but otherwise arbitrarily (see FIG. 
[2]). Then for any state as in (4.16) characterised by its respective Wigner function W(x,p) 



FIG. 2: The maximal support of the function f(x,£) in the coordinates (x, £) for 

b\ = —1, a 2 = 3. 


one has that 




dpW(x,p)A(x,p ) = 0. 


The crucial role of the interference component of a Wigner function will be illustrated in 
the case of the ground state of a 1-D harmonic oscillator. The wave function of this state 

is 


V’o(aj) = 


Moj 


nh ) 


1/4 


exp 


Mux' 

2 h 


2 


and it can be written as the sum ^o(x) = , ipo(-)(x) + ' l l J o {+)( x ) with 


(4.18) 





/ Mux 2 \ 

\~2n~) 


Y(-x) , ^ O (+)(a0 



/ Mux 2 \ 

\~2n~) 


Y{x). 
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On the other hand the Wigner eigenfunction of the ground state of the 1-D harmonic 
oscillator is given by 

1 / p> + MW\ (4ig) 


W 0 (x,p) = — exp 
7m 


KMuj 


which is of the Gaussian type and therefore nonnegative (see 1231 ). 


The figures illustrating: the complete Wigner function (4.19) in the area [0, +oo) x R, the 
Wigner function corresponding to the wave function V’o(+) (#) and the interference Wigner 
function between -0 O (_)(a;) and i/} 0 ^(x) restricted to the region [0, +oo) x R are presented 


at FIGS. 3a, 3mandwc 



(a) The complete Wigner eigen- (b) The Wigner energy eigen- (c) The interference Wigner 
function function without the interfer- eigenfunction 

ence contribution 

FIG. 3: Components of the Wigner eigenfunction in the area [0, +oo) x R for the ground 
state of the 1-D harmonic oscillator. 


For x = 0 the component of the Wigner eigenfunction without the interference part 


disappears. This observation is in agreement with the formula (4.9) applied with the 
parametres ai = 0, bi = oo. 

Although the complete ground state Wigner eigenfunction is positive, its interference 


component admits negative values as can be appreciated at FIG. 3c 


The Wigner function related to the wave function (4.7) represented by a sum of spatially 


separable functions can be also found in another way. For every function 'ipEaibi(x) = 
Y(x — aiY'pEaib, ( X )Y(bi — x) we introduce a new function 

i>Ea l b l (x ) := Y(ai - x) + Y(x - a{)ij} E aibi{x)Y(b t - x) + Y(x - bi), 


which is not square integrable, but it is a tempered distribution. Now the complete energy 
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eigenfunction is the product of functions 

k 

^e{x) = W^EaibXx) , ~oc < ai < bi = a 2 < b 2 = a 3 < ... < b k _ x = a k < b k < oo. 
1=1 


Therefore the Wigner energy eigenfunction can be calculated as the convolution (4.4) 
of Wigner functions representing However, these Wigner functions are usually 

tempered generalised functions and are considered at Appendix [Aj 


4.3. Compatibility conditions for Wigner functions 

In the previous subsection the problem of relations between partial Wigner functions 
defined at spatially separable intervals has been discussed. Now we propose a way in 
which these partial functions can be put together. The lacking elements are coefficients 
standing at them. There are a few ways, in which they can be found. The most natural 
seems to this one based on considerations devoted to an analogous problem in Sec. [3j 
First in each region we solve the equation for phase. For the intervals, at which the series 


expansion (3.8) holds, this is the system (3.9). Close to turning points it is necessary 


to deal with the Schroedinger equation for a polynomial potential. Then from the re¬ 
quirement of continuity for the wave function we obtain factors standing at partial wave 
functions. These factors are transferred to the partial Wigner functions and interference 


Wigner functions with the use of the Weyl correspondence (2.5). 


5. EXAMPLES OF THE WKB CONSTRUCTION IN DEFORMATION QUAN¬ 
TIZATION 

This section contains two examples of application of the WKB method proposed before to 
1-D energy eigenvalue problems in deformation quantization. The first one corresponds 
to the bound states of a 1-D harmonic oscillator and the second to the unbound states of 
a Poeschl - Teller potential. Since these two problems can be solved strictly they provide 
us the possibility to compare exact results with the ones obtained in the first order of the 
WKB approximation. 
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5.1. The WKB approximated Wigner energy eigenstates for the 1—D harmonic 
oscillator 


The Hamilton function of the 1-D harmonic oscillator is given by the formula 

p 2 Mu 2 x 2 


H(x ' p) = 2M + 


(5.1) 


Applying the energy quantization rule (3.28) we obtain that energy levels are determined 
by the relation E n — hu (n + |) , n — 0,1, 2,... exactly like in the strict procedure. The 
corresponding turning points are given by — Xq = — \] an d x ° = \J ■ 

Then combining the applicability condition (3.13) and a demand of compatibility of a 


linear approximation of the potential 


Mu 2 x 2 


with this real potential near to the turn¬ 


ing points we see, that the quasi-classical method can be definitely applied for energy 
eigenvalues fulfilling the condition n > 8. 


Moreover, we divide the domain M. of Eq. (3.7) for the harmonic oscillator into five regions: 

1. the classically forbidden areas A L = (—oo, — |xo] and A R = [|x 0 , oo) in which the 
series expansion (3.8) works, 

2. the intervals B L = [—|x 0 , — |x 0 ] and B R = [|x Q , |(c 0 ] containing the turning points 

3. and the classically accessible part C = [— |xo, |xo] • 

The notation used above is in agreement with the one proposed for FIG. [l] 


In the areas Al,C and A R we solve Eqs. (3.9) up to the first order approximation. In 
the regions B R and B R we deal with the Schroedinger equation with a linear potential. 

A list of phases and normalising factors for the five enumerated regions is presented below, 
where for simplicity we have omitted the index ‘(1)’ at normalising factors. 


1. In the interval A R it is found that 
iMuj 


(TO = 


—x\x* — 


H(2n + 1) h(2n + 1) 


Mu 


Mu 


In 


—x + 


2 _ ft(2rt+l) \ 
x Muj 


h(2n-\-l) 

Muj 




ai — — — In ( k M 2 u 2 x 2 — MJlu(2n + 1)) , 


and Da l = y, where A^ is a parameter determined by the normalisation require¬ 
ment. 
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2. In the region B L the series representation of the phase a does not work. Since the 
potential has been approximated there by a linear function, we see that 


a = —ih In 


- 2 V3^ (2n + 1 ). / Jyi|±t) + 


X 


N 


and the normalising coefficient C Bl = 2 i/ 6 M ^/^/^/\ 2 n+i) 1 / 12 - 

3. In the classically accessible area C there are two physically acceptable solutions. As 
it was shown in Sec. [3j they differ at signs at do but they have the same phases o\. 

Thus 

&oi — —&011 = 

(2 n + 1 )h 


X\ 


Mu 
(2 n + l)h 


Mux 2 
(2 n + 1 )h 


+ arcsm x\ 


Mu 

(2 n + 1 )h 


+ 


(2 n + l)7r 


and 


an — G\ u — — j In ( Mhu(2n + 1) — M 2 u 2 x 2 ) 


The normalising factors are Djc — — and D IIC = 


4. In the region B b one gets 
cr = —ih In I $ 




X 


with the factor C Br = 2 i/e M i /i S/Z^( 2 n+i)^ - 
5. Finally in the interval A B 


iMu 


a 0 = 


x\ x 


h{ 2 n + 1) h{ 2 n + 1) x + 


Mu 


Mu 


In 


2 h(2n-\-l) 


X 


Muj 


h(2n-\-l) 

Mu) 


G\ — — — In ( M 2 u 2 x 2 — Mhu(2n + 1)) 
and the coefficient Da r = . 

Since the explicit formulas for partial and interference Wigner functions obtained via 


the Weyl correspondence (4.1) are extremely complicated, we present the results by the 


following pictures obtained numerically. 
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As it can be seen comparing FIG. 4a and FIG. 4b for n — 8, the first order of the WKB 
approximation applied to the wave energy eigenfunction generates a function which is in 
a good agreement with the strict solution. 



(a) The first WKB approximation (b) The strict wave function 

FIG. 4: The harmonic oscillator wave energy eigenfunctions in the first WKB approxi¬ 
mation and in the strict version for n = 8. 


An analogous situation takes place for the Wigner energy eigenfunction for n = 8. The 
strict Wigner energy eigenfunction and the Wigner energy eigenfunction in the first order 
of the WKB approximation look practically the same. The exact solution, determined by 
the expression 


W hoj{n+1/2 )(x,p) = 


-l) r 


7T h 


exp 


■SW£). 


(5.2) 


where H denotes the Hamilton function (5.1) and L n (y) := ]C!!,=o( — 1) 


n 


n — m 


y 

ml 


is the nth Laguerre polynomial, can be seen at FIG. [5} 
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FIG. 5: The strict Wigner energy eigenfunction for n = 8. 

The effect of the interference terms on the complete Wigner energy eigenfunction can 
be appreciated at FIG. [6j Values of this Wigner interference function are in general 
small comparing to the complete Wigner function but at some points the interference 
contribution is significant. 



FIG. 6: The contribution due to the interference between the wave functions from the 
regions A L and A R to the Wigner function. 

The Wigner energy eigenfunction without any interference contribution is presented at 
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FIG. [FJ This incomplete function is in fact sufficient to deal with the spatial probability. 



-4 


FIG. 7: The Wigner energy eigenfunction without the interference contribution. 

To complete our analysis it is important to mention the following remarks about the Oth 
approximation of the Wigner energy eigenfunction. In this case the phase <Jo is a solution 
of the classical Hamilton - Jacobi equation. However, the class of admissible solutions is 
larger than in the classical physics, because we accept also imaginary phases in classically 
forbidden regions. 

The Wigner energy eigenfunction in the zeroth WKB approximation for the n — 8 is 
presented at FIG. [8} In general it is similar to the strict Wigner energy eigenfunction 
presented at FIG. [5j Since it takes also negative values, it does not represent a probability 
distribution and cannot be treated in a classical way. 
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FIG. 8: The Wigner energy eigenfunction for n = 8 in the zeroth WKB approximation. 


5.2. The unbound energy eigenstates of the Poeschl Teller potential in deforma¬ 
tion quantization 


As the second example we consider a case of unbound states, for which the series expansion 


(|3.8|) is valid everywhere. The chosen potential is the Poeschl - Teller potential described 

h 2 a 2 1 


by the expression 


V{x) = 


M cosh 2 (ax)’ 

where a > 0 is a parameter. The form of this potential is shown at FIG. [9j 


(5.3) 



FIG. 9: The Poeschl - Teller potential for the parametres h — a — M — 1. 
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The energy eigenvalue problem for this potential is solvable for any positive energy E > 0 
and its solution is the function 


, , v , ,'ik - atanh(ax)\ , , 2ME 

ipk{x) = A ( -——-—- exp (ikx) , k = \j ——. (5.4) 


ik + a 


h 


An interesting feature of the potential (5.3) is that every incident particle, regardless of 


its energy, passes right through. The wave function (5.4) is not normalisable and the 


relative spatial density of probability determined by it is illustrated at FIG. 10a 



(a) The spatial density of probability (b) The momenta density of probability 


FIG. 10: The densities of probability for the unbound state for k = 6 and M — a — 1 in 
the Poeschl - Teller potential. 


In the momentum representation the wave function 'ipk(p) is a sum of two generalised 
functions. Indeed, 


ipkip) — A\/2iih- 


ik 


ik + a 


5(p — kh ) 


2 kh 


vp- 


sinh 


r(p—kh) 
2 ah 


(5.5) 


where vp denotes the principal value of the function. The density probability with respect 


to momenta is not defined but the relative probability is illustrated at FIG. 10b 


In the quasi - classical approximation the series expansion (3.8) can be applied everywhere. 
Moreover, the Poeschl - Teller potential is reflectionless and it is assumed that a source 
of particles is localised at — oo. As the limit in the integral determining the phase <7o we 
put xq = 0. Thus 


(Tq = 


Ti\/k 2 cosh 2 ax + 2d 2 
\/k 2 cosh 2 ax + 4a 2 + k 2 


2 arctan 


2 a sinh ax 


a Jk 2 cosh2aa; + 4a 2 + k 2 


+ 
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k ( k sinli ax 

—arcsmh — = 

a \y/2a 2 + k 2 


and 


(5.6a) 


ci — —- In (h cosh ax\Jk 2 cosh 2 ax + 2a 2 'j . (5.6b) 

Since the integral (4.1) is divergent in both: the strict (5.4) and the approximated (5.6a[ 


5.6b) case, to calculate respective Wigner functions is necessary to apply the formula 


(A. 10). Thus to show that the WKB method leads to a reasonable result, we compare 


the wave function (5.4) and the approximated function 2.5-exp (|cr 0 + on) . The values of 


the approximated function have been rescaled by the factor 2.5 because the most natural 
choice of the free term equal to 0 in ctq gives a value of the factor |A| ^ 1 in the formula 


(5.4) 



FIG. 11: Comparison between the real and imaginary parts of the strict wave function 
and its first WKB approximation for the unbound state k = 6, M = a = 1. 


The results can be appreciated at FIGS. |lla| and |llb[ where the strict function is in 
blue and the approximated one in yellow. It is hard to distinguish between them. The 
respective Wigner function can be found at the end of Appendix [B] at FIG. 12 


6. CONCLUSIONS 

In this work a careful analysis of obtaining an approximate Wigner function by means of 
WKB approximation of the wave function is presented. In particular, expressions for the 
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Wigner function coming from a product and a superposition of functions are found. The 
WKB algorithm adapted to the deformation quantization formalism enables us to find 
approximate energy levels and approximated Wigner energy eigenfunctions. Moreover, 
the nth approximation of the Wigner energy eigenfunction is determined by the nth 
approximation of the phase a(x). 

The WKB approximation in deformation quantization is an approach alternative to the 


wave quantum mechanics version of the quasi - classical method. It starts from Eq. (3.5) 


and through formula (4.1) leads directly to the Wigner function so there is no step, at 
which any use of a wave function would be required. 

The initial step of the quasi - classical method is solving of the classical Hamilton - 
Jacobi stationary equation. However it does not mean that in the Oth approximation 
we deal with a ‘classical’ limit of the Wigner function. The reason is that in quantum 
considerations also imaginary solutions of the Hamilton - Jacobi equation are acceptable. 
There are some important difficulties when the quasi - classical method is applied. The 


principal one arises from the non - locality of the Weyl correspondence (2.5). Indeed, 
to find the Wigner energy eigenfunction at an arbitrary point (x,p) of the phase space 
we have to know the complete phase (phases) a(x) or, equivalently, the complete wave 
function. Now, since the wave function is in general a sum of spatially separated functions, 
interference terms are present and for a problem with several spatial regions it is time - 
consuming to calculate the complete Wigner function. 

Another disadvantage of the proposed algorithm is the fact, that on the contrary to the 
idea of the Hamilton formalism, the considerations have been done from a spatial point 
of view. Thus positions and momenta are not treated on an equal footing. 
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Appendix A: The Wigner function of a nonnormalisable state 


In the case when a wave function is not square integrable, the integral in formula (2.5) of 


a Wigner function may not be convergent. Thus it is necessary to propose some extension 


of expression (2.5) applicable to states represented by generalised functions. 


To establish notation we quote three definitions: 

Definition A.l. Let ip(z) be a function from the Schwartz space S. Then its inverse 
Fourier transform is determined by the integral 

1 


FlvizM) = Fz[<p](t) ■■= 


<p(z) exp (—izt)dz. 


a/2vt J -o 

Generalisation of Def. |A.1| on functions of many variables is straightforward. 


(A.l) 


The inverse Fourier transform of a function from the Schwartz space is also an element of 
S. A tempered distribution T is a linear continuous functional over the vector space S. 
The set of tempered generalised functions will be denoted by S'. 

Definition A.2. An inverse Fourier transform J~[T] of a tempered distribution T is a 
tempered generalised function satisfying the equality 




(A. 2) 


for every <p 6 S. 


Definition A.3. The convolution of two tempered generalised functions S and T (if it 
exists) is a generalised function S *T defined by the condition 


V cp e S (^S *T,ip^ = (^S x <8)T y ,(p(x+ y)Y 


(A.3) 


The inverse Fourier transform of the product of tempered distributions S ■ T (if it is well 
defined) is the convolution of their respective inverse Fourier transforms 


T[S-T]=-j=F[S]*m. 


(A.4) 


It is well known that the inverse Fourier transform satisfies the following properties: 
1. a translation of the argument leads to the relation 


V a G M iF[T(z + a)](t) = exp (iat)T[T(z)](t), 


(A.5a) 
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2. a multiplication of the argument by a real number implies the equality 


v a € R \ {0} nT(az)](t) = A./[T( Z )] ( - ) , 

a V a, 


(A.5b) 


3. the complex conjugation of a tempered generalised function T is transformed into 


F{T(z)](t)=Tp ’(*)](-«), 

4. the inverse Fourier transform of the derivative 

dT(z 


(A.5c) 


T 


dz 


(t) = itr\T(z) ](«). 


(A.5d) 


One sees that Def. (2.5) of a Wigner function related to a state |"0) can be written as 

Te — ~) (v), (A.6) 


W(x,p ) = 




P\ 

h) ’ 


where £ is a variable and x plays a role of a parametre. We do not assume that |-0) is an 
energy eigenstate. 

Applying the fact that the inverse Fourier transform of the product of functions is the 
convolution of inverse Fourier transforms one obtains that 


W(x,p ) = 


2nh 




i/j (x + ^ 


* J-e 


( x- | 


(A.7) 


Using properties (A.5a), (A.5b) and (A.5c) one gets 


V’ ( x + | 


{t) = 2 exp (2 itx) -FfoKO] ( _2t ) 


and 




( x- | 


(t) = 2 exp (-2 itx) (-2*0 • 

Therefore the Wigner function representing a state | -0) equals 

Vh(x,p) = —r: fexp(ita).F[^(f)](£) * t exp(ita).F[^(f)](*)) (-y- 


or 


W(x,p) = ^ + z)](t) * t J r [^(^ + x)](t)j I ■-~ 


(A.8a) 


(A.8b) 


(A.9) 


(A.10) 
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Appendix B: The Wigner function of an unbound state in the Poeschl —Teller po¬ 
tential 


From the relationship (A.9) one can see that the starting point for calculating a Wigner 
function for a state 1-0) is the inverse Fourier transform of the wave function of this 


state. In our problem the wave function is expressed as (5.4). Direct calculations lead to 
the conclusion that 


•Fx 


. , ik — atanh(ax) , . , . 

A ( ---- —exp (ikx) 

ik + a 


(t) = ——— \j— ( 2k5(t — k) — vp 


a + ik\ 2 


sinli 


o-b 

2 a 


(B.l) 


Therefore the Wigner eigenfunction for an eigenvalue E = FF- equals to the following 


convolution 


W h 2 k 2 (x,p) = 




2 h(a 2 + k 2 ) 


2keM-ikxW-k)-vp eXp( ~ te) N |», 


2k exp(ikx)S(t — k) — vp- 


exp (itx) 


sinh 


r (k—t) 
2 a 


sinh 

2 p 

h 


t (k—t) 
2 a 


(B.2) 


The most complicated part is to derive the convolution 


exp (—itx) exp (itx) 

vp-—-r- *t vp 


sinh 


r (k—t) 
2 a 


sinh 


( fc -b 

2 a 


(z). 


The linear change of the variable k — t = u turns the above formula into 

(2k — z). 


exp (ixu) exp (—ixu) 

V P . , Wvp- 


sinh (=) sinh (g) 


(B.3) 


Its imaginary part vanishes so finally the expression (JB.3|) consists of two components 
cos (xu) 


COS OT 


vp 


sinh (!) VP sinh (f) 


(2k — z) + 


sin ot 


sin xm 


s in h m U sinh (f) 


(2k — z). 


The second element is a convolution of functions and can be found by integration as 

f + °° sin(xu) sin(x[2/c — z — u]) 


sinh (^ 


( 2 a ) sinh ^ 


n[2k—z—u] 
2 a 


-du. 


(B.4) 
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On the other hand 


exp(ix«) exp (ixu) 

Re < vp- t— r- * u vp- t— v 

' sinh(g) P sinh(g) 

Therefore 

exp (ixu) exp {—ixu) 


cost xu) cos [xu) sin(xw) sin(xM) 

= V P ■ , *nW — - - - - .-x *> 


S mh(g) ‘sinh(g) sinh(g) “sinh(g)' 


vp 


sinh S) *“ vp sinh ®) 


„ . exp (ixu) exp (ixu) 

= Re { V P ■ , / tru \ *» V P 


sinfxitj sm(OT) 

+2 :—7 


sinh(^) " sinh(f)J sinh (f) “sinh (f^) 


But 


Indeed, 


exp (hen) 
VP sinh (f) 

r ‘ +0 ° cos tu 




l~2 

iay — tanh (a(f + x)) 


(«)■ 


r»+oo 


vp 


^ (g) 


du = 0 and 


sin tu 

™h(g) 


du = a tanh at 


see im The observation ( |B.5[ ) implies that 

exp (ixu) exp (ixu) 9 j2 ~ r l9/ . , N1 . . 

vp ^h7=T *“ vp ^hM = “ 2o V ^*[ tanh (o(i + x)) ] (M) ' 


:g) ^(S) 

Moreover, the derivative 

d tanh (a(t + x)) 
dt 


= a (l — tanh 2 (a(t + x))) . 


Applying property (A.5d) to the expression (B.6) we obtain that 


exp(ixn) exp(zxn) 
VP sinh (g) * u VP sinh (f) 


(z) = -4 a 2 5(z) + 2z • vp- 


exp(ixz) 


sinh (i) 


Thus 


exp(ixn) exp (—ixu) 
VP sinh (g) *“ VP S inh(g) 


{z) = —Aa 2 S(z) + 2cos(xx) 


sinh (i) 


sin(xn) sin(xn) 


si„h(g) u sinh (g) 


(7). 


Finally 


, 2 k 2 - a 2 2\A\ 2 k ( kh + p 

W h 2 k 2 (x,p) = \A\ 2 — —— -5(p + hk) - cos ( 2x 


k 2 + a 2 

2|A| 2 


h{a 2 + k 2 ) 


ft 


vp 


(B.6) 


(B.6) 


+ 


sinh ( 


cos ( 2x 


kh + p\ kh + p 


h 2 {a 2 + k 2 ) V. ft J sinh ^ A^+p) j 
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H 1 2 3 4 5 


L 


+oo 



x{2kh-\-2p—u) 
h 


du. 


(B-7) 


h(a 2 + B) 


—oo 


( 2 a ) sinli ^ 


7r(2kh-\-2p—u) 
2ah 


The unusual relation that the momentum p corresponds to the wave vector —k results 
from the fact that our sign convention in the * - product is in agreement with Fedosov’s 
papers [211 22] ■ 

Notice that for the state satisfying the condition \k\ — a the component containing the 
Dirac delta disappears. 


-8-7 P 



4 


FIG. 12: The strict Wigner function of the Poeschl - Teller potential for k = 6, M = a = 

1 . 
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